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The  V5+  electrolyte  solution  from  Vanadium  Redox  Flow  Batteries  was  studied  by  variable  temperature 
170  and  51 V  Nuclear  Magnetic  Resonance  (NMR)  spectroscopy  and  density  functional  theory  (DFT)  based 
computational  modeling.  It  was  found  that  the  V5+  species  exist  as  hydrated  penta  co-ordinated  vanadate 
ion,  i.e.  [V02(H20)3]1+.  This  hydrated  structure  is  not  stable  at  elevated  temperature  and  change  into 
neutral  H3VO4  molecule  via  a  deprotonation  process  and  subsequently  leading  to  the  observed  V2O5 
precipitation  in  V5+  electrolyte  solutions. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Recently,  Vanadium  Redox  Flow  Batteries  (VRFB)  have  shown 
great  promise  as  an  efficient  large  scale  energy  storage  system 
for  a  wide  range  of  applications  including  uninterruptible  power 
supply,  stabilization  of  wind  turbine  output,  and  hybrid  photo¬ 
voltaic  cells  [1-4].  Unlike  other  popular  electrochemical  energy 
storage  systems  (such  as  lithium  ion  batteries),  VRFB  store  energy 
in  electrolyte  solutions  rather  than  in  electrodes.  This  setup  could 
result  in  a  considerable  increase  in  life  time  due  to  low  self 
discharge  of  the  battery  provided  that  the  charged  electrolyte 
solutions  can  be  kept  in  stable  conditions.  In  reality,  the  fully 
charged  V5+  electrolyte  solution  displays  poor  stability  at  elevated 
temperatures  (>310K)  and  also  at  high  vanadium  concentrations 
(>2M)  [5-7].  This  poor  stability  is  witnessed  as  hydrated  V205 
precipitation,  which  leads  to  energy  loss  and  failure  of  the  bat¬ 
tery  [8,9].  The  incapacity  to  increase  the  vanadium  concentration 
(>2  M)  in  the  electrolyte  solution  without  increasing  the  sulfuric 
acid  concentration  (>5M)  puts  a  limit  on  the  maximum  energy 
density  of  the  VRFB  (<25  Whkg-1).  Up  to  now,  the  chemistry 
behind  the  precipitation  behavior  of  V5+  electrolyte  solution  is  not 
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well  understood.  Understanding  the  chemistry  behind  this  pre¬ 
cipitation  behavior  is  imperative  for  developing  a  stabilizer  or 
even  a  new  stable  V5+  electrolyte  solution,  which  can  ultimately 
increase  the  life  cycle,  energy  density  and  reduce  the  operating 
cost  of  the  VRFB.  In  this  work,  a  combination  of  Nuclear  Mag¬ 
netic  Resonance  (NMR)  and  density  functional  theoretical  (DFT) 
studies  was  carried  out  to  provide  a  clear  view  about  structure  of 
V5+  species  and  their  precipitation  behavior  at  elevated  tempera¬ 
tures. 

2.  Experimental 

The  V5+  electrolyte  solution  is  prepared  by  electrochemically 
charging  the  electrolyte  stock  solution  (2  M  V0S04xH20  dissolved 
in  5M  sulfuric  acid)  using  a  flow  test  cell  with  the  current  den¬ 
sity  of  1  A  and  voltage  of  1.7  V.  The  5  ml  of  V5+  electrolyte  solutions 
is  taken  from  the  positive  half  cell  after  12  cycles  in  which  10  |xl 
of  commercial  H2170  solution  (20%  enrichment)  was  added  for 
170  NMR  measurements.  The  variable  temperature  170  and  51V 
NMR  measurements  were  performed  using  a  Varian  500  Inova 
spectrometer  (Bo  =  11.1  T  and  170  and  51 V  Larmor  frequencies  are 
67.8  and  130.7  MHz,  respectively)  and  their  chemical  shifts  were 
externally  referenced  to  a  commercial  H2170  and  clear  V0C13  solu¬ 
tions,  respectively.  DFT  calculations  were  carried  out  using  the  ADF 
2009.01  program  [10].  The  hybrid  B3LYP  function  with  QZ4P  ( quad 
Z,  4  polarization  functions,  all  electron)  basis  set  with  the  zeroth 
order  regular  approximation  (ZORA)  is  used  for  both  geometry  opti¬ 
mization  and  NMR  calculation.  All  the  DFT  calculations  reported 
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Fig.  1.  Variable  temperature  chemical  shift  analysis  of  170  NMR  spectra  (left)  and  51 V  NMR  spectra  (right)  of  2  M  V5+  electrolyte  solution  measured  at  a  magnetic  field  of 
11.7T.  The  open  circles  and  the  closed  circles  in  the  chemical  shift  analysis  represent  heating  cycle  and  subsequent  cooling  cycle,  respectively.  The  error  in  chemical  shift 
measurements  is  ±0.5  ppm. 


here  are  carried  out  under  conductor-like  screening  model  (COSMO) 
with  water  as  solvent  [10]. 

3.  Results  and  discussions 

Fig.  1  shows  the  temperature  dependent  isotropic  chemical  shift 
of  170  and  51 V  NMR  spectra  measured  on  freshly  prepared  2  M  V5+ 
electrolyte  solution.  The  chemical  shift  of  both  170  and  51 V  NMR 
shows  an  irreversible  change  around  335  K,  which  indicates  a  ther¬ 
mally  induced  irreversible  structural  transition  of  the  vanadium 
species  in  the  V5+  electrolyte  solution.  This  structural  transition  due 
to  heat  treatment  is  evidenced  by  visible  yellowish  red  particles  in 
the  electrolyte  solution  few  hours  after  the  NMR  measurements 
and  the  precipitated  particles  are  known  to  be  hydrated  V2O5  [11]. 
To  understand  the  pathway  of  V205  precipitation,  the  structures 
of  the  vanadium  species  before  and  after  the  structural  transition 
must  be  identified.  It  is  widely  believed  that  oxovanadium  (V02+) 
cation  is  the  major  species  in  highly  acidic  solutions  (pH  <  2)  such  as 
our  VRFB  electrolytes  [12,13].  This  V02+  species  reportedly  exists 
in  its  hydrated  form  with  a  six-coordinated  octahedral  structure 
in  which  two  oxygen  atoms  are  coordinated  in  cis-configuration 
along  with  four  complexed  water  molecules  (i.e.  [V02(H20)4]+) 
[14].  However,  our  DFT  based  geometry  optimization  revealed  that 
this  octahedral  vanadate  ion  structure  is  not  energetically  favorable 
and  a  slightly  more  stable  structure  is  formed  by  repelling  a  water 
molecule  located  at  trans  position  to  the  vanadyl  oxygen. 

This  turns  the  structure  into  a  penta  co-ordinated  [V02(H20)3]+ 
molecule  (hereafter  called  compound  A)  as  shown  in  Fig.  2,  which  is 
in  good  agreement  with  results  of  Buhl  and  Parrinello  [15].  Hence, 
compound  A  is  the  plausible  initial  structure  in  the  V5+  electrolyte 
solution  before  the  thermally  induced  structural  transition.  Since 
the  structural  transition  is  followed  by  V205  precipitation,  the  most 
probable  way  this  process  can  occur  is  by  condensation  of  mononu¬ 
clear  vanadium  species  containing  hydroxyl  groups  [16].  Therefore 
the  new  vanadium  species  formed  at  high  temperature  (>335  K) 
most  likely  contain  some  hydroxyl  groups  in  its  structure.  The  pos¬ 
sible  scenario  to  form  any  such  hydroxyl  group  is  by  deprotonation 
of  water  molecules  complexed  in  the  initial  compound  A  which 
could  lead  to  a  neutral  species  H3VO4  (hereafter  called  compound 
B)  as  shown  below: 

[V02(H20)3]+^H3V04  +  H30  (1) 


It  is  known  that  such  deprotonation  can  happen  if  there  is  a 
change  in  pH  and/or  a  change  in  the  temperature  of  the  vana¬ 
dium  solutions.  This  neutral  molecule  H3VO4  (or  more  precisely 
VO(OH)3)  is  well  known  to  exist  at  room  temperature  and  at  low 
pH  (~3.5)  [14,17].  Recently,  Sadoc  et  al.  [18]  reported  that  the 
hydrated  V02+  molecule  is  prone  to  deprotonation  at  high  tem¬ 
peratures  (500 1<),  resulting  in  the  formation  of  V0(0H)2(H20)2+ 
and  subsequently  the  neutral  H3V04.  The  optimized  structure  of 
neutral  compound  B  using  ADF  is  shown  in  Fig.  2.  We  hypothesize 
that  deprotonation  in  highly  acidic  medium  (<1  pH),  such  as  our 
V5+  electrolyte  solution,  is  possible. 

To  confirm  our  hypothesis,  we  calculated  the  chemical  reaction 
energy  (A  E  =  EA  -  EB  -  £h3o)  of  proposed  deprotonation  reaction 
(Eq.  (1 ))  using  total  energy  of  DFT  optimized  compounds  A  (EA)  and 
B  (£b)-  The  calculated  A  E  of  Eq.  ( 1 )  shows  it  as  an  endothermic  reac¬ 
tion  (AE~1.2eV).  This  means,  the  proposed  structural  transition 
and  subsequent  precipitation  need  thermal  energy  which  agrees 
well  with  the  reported  poor  stability  of  V5+  of  electrolyte  solution 
at  high  temperatures  [5-7].  Since  this  deprotonation  reaction  is 
endothermic,  it  is  more  likely  to  be  a  reversible  reaction  such  that 
the  resultant  compound  B  can  be  reversed  to  compound  A  by  sim¬ 
ple  protonation  at  any  given  temperature.  The  isotropic  chemical 
shift  (5)  parameters  for  the  isolated  compounds  A  and  B  were  cal¬ 
culated  using  ADF  program  and  compared  with  experimental  NMR 
chemical  shifts.  It  should  be  noted  that  experimental  NMR  chem¬ 
ical  shifts  observed  during  cooling  cycle  represents  the  averaged 
state  of  both  compounds  A  and  B.  Hence  more  attention  should  be 
paid  to  the  trends  in  the  DFT  calculated  8  values  between  the  com¬ 
pounds  rather  than  to  their  actual  values  [15].  The  DFT  computed 
51 V  chemical  shifts  for  compound  A  (5  ~  -623  ppm)  and  compound 
B  (5  ~  -577  ppm)  show  that  any  transition  from  compounds  A  to 
B  should  induce  a  shift  towards  higher  frequency  (i.e.  nuclear  de¬ 
shielding)  in  the  51 V  NMR  spectra,  which  is  in  good  agreement  with 
our  experimental  data  (see  Fig.  1).  Similarly,  the  computed  170 
chemical  shift  for  water  molecules  in  compound  A  {8  =  -23  ppm) 
and  hydroxyl  groups  in  compound  B  (8  =  -410  ppm)  reveals  that 
transition  between  compounds  A  and  B  would  lead  to  significant 
170  chemical  shift  towards  lower  frequency  (i.e.  increase  in  nuclear 
shielding),  which  also  coincides  with  our  experimental  data. 

The  structural  transition  in  Eq.  (1)  is  mainly  driven  by  proton 
migration  between  the  water  and  hydroxyl  molecules  in  com¬ 
pounds  A  and  B.  Such  a  proton  migration  in  the  solution  is 
bound  to  directly  influence  the  170  NMR  line  width  [19,20].  Fig.  3 
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Fig.  2.  Geometry-optimized  structures  for  [V02(H20)3  ]+H20  (compound  A)  and  the  neutral  VO(OH)3  (compound  B).  The  vanadium,  oxygen  and  proton  atoms  are  represented 
as  pink,  red  and  white  spheres,  respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  caption,  the  reader  is  referred  to  the  web  version  of  the  article.) 


shows  170  line  width  during  heating  cycle  (only  compound  A) 
and  cooling  cycle  (representing  both  compounds  A  and  B).  The 
relatively  smaller  line  width  during  cooling  cycle  indicates  the 
possible  faster  proton  exchange  in  solution  and  if  this  fast  proton 
exchange  happens  directly  between  hydroxyl  groups  of  two  com¬ 
pound  B  molecules,  it  will  facilitate  the  condensation  followed  by  a 
hydrolysis  reaction  (2H3V04  -*  V205  +3H20)  [16].  Such  a  hydroly¬ 
sis  reaction  can  lead  to  the  formation  of  V-O-V  bonding  groups 
known  as  structural  building  mosaics  of  the  solid  V2O5  which 
would  ultimately  lead  to  the  crystallization  of  V205  [16].  Further, 
it  is  interesting  to  note  that  the  170  NMR  spectra  shows  rela¬ 
tively  larger  line  width  at  low  temperatures  (<0  °C)  which  indicates 
slower  proton  exchange  rates  in  the  solution  at  this  low  tempera¬ 
ture  region.  This  observed  slow  proton  exchange  proves  that  there 
is  no  significant  protonation/deprotonation  reaction  at  low  temper¬ 
atures  and  explains  the  relatively  better  stability  of  V5+  electrolyte 
solutions  at  low  temperatures. 


Temperature  (K) 


Fig.  3.  Variable  temperature  170  NMR  line  width  analysis  of  2  M  V5+  electrolyte 
solution  measured  at  11.7T.  The  open  circle  and  the  closed  circle  represent  heat¬ 
ing  cycles  and  the  subsequent  cooling  cycles,  respectively.  The  error  in  line  width 
measurements  is  ±100  Hz. 


Now  the  question  arises,  why  the  other  vanadium  species  (i.e. 
Vz+,  z  =  2-4)  despite  their  similar  hydrated  structures  as  compound 
A  have  more  stability  and  show  no  structural  transition  associated 
with  deprotonation  [21].  In  the  electrolyte  medium  the  vanadium 
cations  are  solvated  by  water  molecules  forming  hydrated  cations 
[V(H20)n]y+  where  the  electron  density  from  molecular  orbital  of 
the  water  molecules  is  transferred  to  the  empty  orbitals  of  the 
metal  cation  [16].  This  charge  transfer  weakens  the  O-H  bond  in 
the  coordinated  water  molecules  and  makes  it  more  acidic  leading 
to  deprotonation  of  the  hydrated  vanadium  cations  as  follows  [16]: 

[Vz+(H20)n]y+  +  /lH20  [Vz+(0H)f|(H20)n_f,](y-h)+  +  /1H3O+  (2) 

In  this  deprotonation  reaction,  the  hydrolysis  ratio  (h),  i.e.  the 
number  of  hydroxyl  groups  created,  is  directly  proportional  to  the 
charge  of  the  vanadium  cation  (z+)  and  inversely  proportional  to 
the  pH  of  the  solution,  which  is  explained  in  the  charge-pH  dia¬ 
gram  reported  by  Jorgensen  [16,22].  In  the  case  of  compound  A, 
the  higher  charge  of  the  vanadium  cation  (V5+)  in  combination 
with  the  increase  in  temperature  leads  to  a  maximum  hydrolysis 
ratio  (h  =  3 )  leading  to  the  neutral  compound  B.  For  other  vanadium 
species  (z+<  5)  the  lower  charge  prevents  the  total  hydrolysis  of 
the  vanadium  cation  and  provides  greater  stability  of  the  hydrated 
structure.  It  is  reported  that  the  increase  in  the  sulfuric  acid  con¬ 
centration  in  the  electrolyte  medium  increases  the  stability  of  V5+ 
species  [6].  This  can  be  explained  based  on  the  above  discussion, 
as  the  increase  of  sulfuric  acid  reduces  the  total  pH  of  the  solu¬ 
tion  and  hence  decreases  the  deprotonation,  thereby  stabilizing  the 
V5+  species.  This  observation  coincides  with  above  discussion  and 
provides  additional  proof  for  our  hypothesis. 

4.  Conclusions 

The  V5+  cation  structure  and  its  thermal  stability  in  electrolyte 
solution  are  studied  using  170  NMR,  51V  NMR,  and  molecular 
modeling  techniques.  The  V5+  species  exist  as  hydrated  penta  co¬ 
ordinated  vanadate  ion,  i.e.  [V02(H20)3]1+  in  electrolyte  solution. 
This  hydrated  structure  is  found  to  be  stable  in  the  low  tempera¬ 
tures  (<330  K)  but  at  elevated  temperature  it  changes  to  neutral 
H3VO4  molecule  by  deprotonation  process.  This  neutral  species 
undergoes  a  condensation  reaction  through  its  highly  reactive 
hydroxyl  groups  and  forms  V2Os  precipitate  in  the  electrolyte  solu¬ 
tion.  Hence,  the  precipitation  can  be  avoided  by  preventing  the 
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deprotonation  process  which  is  followed  by  the  structural  tran¬ 
sition  to  neutral  H3V04  molecule.  A  simple  way  to  prevent  the 
deprotonation  process  is  to  increase  the  solvent  acid  concentra¬ 
tion  and  this  explains  the  previously  reported  stabilization  of  V5+ 
electrolyte  solution  at  high  acid  concentrations.  Overall,  it  is  neces¬ 
sary  to  stabilize  the  hydrated  vanadium  cation  structure  to  achieve 
greater  stability  of  the  V5+  electrolyte  solution. 
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